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Abstract 

In this paper we use a generalization of the Voronoi partition to formulate and solve a heterogeneous distributed locational 
' optimization problem for autonomous agents having limited range sensors. Agents equipped with sensors having heterogeneity 
' in their capabilities, communication equipment, and computational capability are to be optimally deployed in a domain of 
\ interest. The optimal deployment is found to be a variation of the generalized centroidal Voronoi configuration, where the 
, sensors are located at the centroids of the corresponding generalized Voronoi cells. We provide a few formal results on stability, 
• convergence, and on spatial distributedness of the proposed control laws under some constraints on the agents' speeds such as 
' limit on maximum speed and constant speed. We support the theoretical results with illustrative simulation results. 
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• 1 Introduction 

1.1 Multi-Agent Systems 

] Technological advances in areas such as wireless com- 
. munication, autonomous vehicular technology, compu- 
' tation, and sensors, faciUtate the use of large number 
I of agents (UAVs, mobile robots, autonomous vehicles 
■ etc.) , equipped with sensors, communication equipment, 
and computation ability, to cooperatively achieve var- 
ious tasks in a distributed manner. Distributed multi- 
agent systems have been shown to achieve and main- 
tain formations, move as flocks while avoiding obstacles, 
etc., thus mimicking their biological counterparts. They 
can also be used in applications such as search and res- 
cue, surveillance, multiple source identification, and co- 
operative transportation. The major advantages of dis- 
tributed systems are immunity to failure of individual 
agents, their versatility in accomplishing multiple tasks, 
simplicity of agents' hardware, and requirement of only 
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minimal local information. At the same time it is impor- 
tant to design distributed control laws that guarantee 
stability and convergence to the desired collective be- 
havior under limited information and evolving network 
configurations. One very useful application of the multi- 
agent systems is sensor network, where a group of au- 
tonomous agents perform cooperative sensing of a large 
geographical area. In this paper, we address the problem 
of optimal deployment of autonomous agents equipped 
with sensors, communication equipment, and computa- 
tional capability. 

1.2 Related Literature 

Advances in the fields of wireless communication, sen- 
sors, computation etc., have led to increased research 
interest in the area of sensor networks, where a large 
number of sensors with limited communication capabil- 
ities are deployed in the domain of interest. Cassandras 
and Li [1] provide a survey on sensor networks. An ex- 
ample of such a sensor network is a network of satel- 
lites equipped with imaging sensors and used to obtain 
the map of a large geographical area where each satel- 
lite provides the map of a small area. Li and Cassan- 
dras [2] represent frequency of occurrence of an event as 



Preprint submitted to Automatica 



25 August 2009 



a density function. The network model is probabilistic 
and assumes that the sensors make observations inde- 
pendently and maximize the expected event detection 
frequency, incorporating communication costs. Zou and 
Chakrabarty [3] use the concept of virtual force to solve 
a similar problem. Hussion and Stipanovic [4]- [5] define 
the effectiveness of coverage so as to ensure that at least 
C* measurements are made at each point in the mission 
domain, where C* is a specified threshold on a suitably 
defined coverage measure. In [5], they ensure collision 
avoidance, while flocking behavior with collision avoid- 
ance is addressed in [6]. Hussein [7] uses Kalman filter 
for multi-agent coverage, in a setting where the sensors 
have noise. 

One of the main problems addressed in sensor networks 
is optimally locating the sensors so as to maximize the 
measurement quality. This class of problem belongs to 
the problem of locational optimization or facility loca- 
tion [8,9]. A centroidal Voronoi configuration is a stan- 
dard solution for this class of problems [10] . Voronoi de- 
composition or Dirichlct tcsselation is a widely used par- 
titioning scheme. It finds application in image process- 
ing, CAD, sensor coverage, multi-agent search [11] and 
many more areas. Cortes et al. [12] address the prob- 
lem of optimal deployment of sensors with limited range 
in a spatially distributed manner using Voronoi parti- 
tion. Pimenta et al. [13] follow an approach similar to 
[12] to address a problem with heterogeneous robots. 
The authors let the sensors be of different types (in the 
sense of having different footprints) and relax the point 
robots assumptions. Power diagram (or Voronoi diagram 
in Laguerre geometry), a generalization of the standard 
Voronoi partition, is used to account for different foot- 
prints of the sensors (assumed to be discs). Due to as- 
sumption of finite size of robots, the robots are assumed 
to be discs and a free Voronoi region that excludes the 
robot footprints, is defined and a constrained locational 
optimization problem is solved. These authors also ex- 
tend the results to non-convex environments. In [15], 
authors consider agents with sensors having anisotropic 
effectiveness, in the sense that effectiveness of the sen- 
sor at a point in space depends on the direction along 
with the Euclidian distance from the sensor. These au- 
thors use a non-Euclidian distance measure and use an 
anisotropic Voronoi partition to solve the problem of op- 
timal sensor deployment. 

1.3 Main contributions 

Most authors consider sensor network to be homoge- 
neous in nature. Whereas, in practical problems, the sen- 
sors may have different capabilities even though they are 
similar in their functionality. The heterogeneity in capa- 
bilities could be due to various reasons, the chief being 
the difference in specified performance. In this paper we 
address a locational optimization problem for sensors 
with heterogeneous capabilities. 



Voronoi decomposition is one of the tools used in loca- 
tional optimization problems. But the existing Voronoi 
decomposition scheme and its variations cannot be used 
for solving the heterogeneous locational optimization 
problem. In this paper, we propose a generalization of 
Voronoi partition, based on the standard Voronoi par- 
tition and its variations. Here a concept of node func- 
tions is introduced in place of the usual distance mea- 
sure. We use this generalized Voronoi partition to for- 
mulate the heterogeneous locational optimization prob- 
lem. The mobile sensors are assumed to have heteroge- 
neous capabilities in terms of the sensor effectiveness. A 
density distribution is used as a measure of probability 
of occurrence of an event of interest. We show that the 
optimal deployment configuration is a variation of the 
centroidal Voronoi configuration. We propose a propor- 
tional control law to make the sensors move toward the 
optimal configuration. Assuming first order dynamics 
for the mobile sensors, we prove, using LaSalle's invari- 
ance principle, that the trajectories of the sensors con- 
verge to the optimal configuration globally and asymp- 
totically. We further analyze the problem in presence of 
some constraints on the agents' speeds and with limit on 
sensor range. Some preliminary results on heterogeneous 
limited range locational optimization problem were re- 
ported in [14] 

1.4 Organization of the paper 

In Section 2, we provide a few mathematical concepts 
used in this paper. We propose a generalization of 
Voronoi partition in Section 3. Section 4 introduces 
the distributed heterogeneous locational optimization 
problem along with the objective function, its critical 
points, the control law and its properties. In Section 5 
we analyze the problem with limit on agents' maximum 
speed and with a constraint on agents to move with a 
constant speed initially and slow down as they reach 
the critical points. We consider the sensor range limits 
and address the limited range distributed heterogeneous 
locational optimization problem in Section 6. Section 7 
provides illustrative simulation results and the paper is 
concluded in Section 8 with some discussions. 

2 Mathematical preliminaries 

In this section we preview a few mathematical con- 
cepts from graph theory, spatially distributed functions, 
LaSalle's invariance principle and Liebniz theorem used 
in the present work. 

2. 1 Concepts from graph theory 

Graphs are extensively used in analysis of multi-agent 
systems and sensor networks. Graph theory [20] provides 
an excellent tool to represent connectivity of agents or 
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sensors in a network. In this section we preview a few 
important terminologies used in graph theory. 

A graph Q = (14, £) consists of a vertex set U, and an 
edge set f G 2^^^ . A graph Q is said to be undirected if 
(i, j) G f ^ (j, i) G E. The map Mg :U ^ 2^ associates 
the set Ng (i) of its neighbors in Q with vertex i. A graph 
is complete if £ = U xU\ diag{U x U), where diagfJA x 
U) = {(u, w) € U xU} . A path connecting vertex i to j is 
a sequence of vertices {iq = . . . ,ik, ik+i — j} with 
the property that G £ for all ? G {0, . . . , A:}. A 

graph G is said to be connected if there exists a path 
connecting any two vertices. 

2.2 Spatially- distributed functions 

One of the desirable properties of a multi-agent system 
is spatial distributedness. A centralized system requires 
that all agents are accessible to the central controller 
and the failure of the central controller leads to failure 
of the entire system. Here, we formally define spatially- 
distributed functions. 

A function / : 1— > is spatially distributed 
over graph Q, if there exists fi : Q x P 1— > K, for 
« G {1,2, . . . ,iV}, such that 



(iv) N is the largest invariant set in E. 

Then, every solution of (1) starting from a point in M 
approaches N as t ^ 00. 

Here, by invariant set we mean that if the trajectory is 
within the set at some time, then it remains within the 
set for all time. Important differences of the LaSalle's 
invariance principle as compared to the Lyapunov theory 
are (i) V is required to be negative semi-definite rather 
than negative definite and (ii) the function V need not 
be positive definite (see Remark on Theorem 3.8 in [23] , 
pp 90-91). 

2.4 Leibniz theorem and its generalization 

The Leibniz theorem is widely used in fluid mechanics 

[26], and shows how to differentiate an integral whose 
integrand as well as the limits of integration are functions 
of the variable with respect to which difiierentiation is 
done. The theorem gives the formula 
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fm = h{puMg{Pi,V)) 

for every V G P„ . Here fi is the i-th component of / and 
Y is the range of fi. 

For a spatially distributed function, the local informa- 
tion is suSicient to compute its value at any given node 
in the graph. That is, the information available from the 

given node itself, and that from its neighbors Afg{pi, T), 
is sufficient to evaluate the function at that node. 

2.3 LaSalle's invariance principle 

Here we state LaSalle's invariance principle [21,22] used 
widely to study the stability of nonlinear dynamical sys- 
tems. We state the theorem as in [23] (Theorem 3.8 in 
[23]). 

Consider a dynamical system in a domain D 

X = fix), f-.D^R" (1) 

Let V : D ^ Rhe a, continuously differentiable function 
and assume that 

(i) M C -D is a compact set, invariant with respect to 
the solutions of (1). 

(ii) y < in M. 

(in) E : {x : X G M, and V{x) = 0}; that is, E is set of 
all points of M such that V{x) =0. 



Eqn. (2) can be generalized for a d-dimensional Eu- 
clidean space as 

4- I F{x,y)dV= f ^dV+ I n{x).u{x)FdS 

dy Jv{y) Jv dy Js 

where, V C R is the volume in which the integration is 
carried out, dV is the differential volume element, S is 
the boimding hypersurface ofV, n(a;) is the imit outward 
normal to S and u(a;) = ^{x) is the rate at which the 
surface moves with respect to ?/ at a; G S. 



3 Generalization of the Voronoi partition 

Voronoi partition [16,17] is a widely used scheme of par- 
titioning a given space and finds applications in many 
fields such as CAD, image processing and sensor cover- 
age. We briefly preview the standard Voronoi partition 
and then present a generalized Voronoi partition used in 
this work. 

3. 1 Standard Voronoi partition 

A collection {Wi}, i G {1,2,..., A''} of subsets of a space 
X with disjoint interiors is said to be a partition of X if 

yJ^W^ = X. 

Let Q C R'', be a convex polytope in d-dimension Eu- 
clidean space. Let V = {pi,f>2, • • • ,Pn}, Pi G Q, be the 
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set of nodes or generators in Q. The Voronoi partition 
generated by V with respect to an Euchdean norm is the 
collection {Vi('P)}jg{i 2,...,^}; and is defined as, 

Vi{V) = {q&Q\ \\q-p,\\<\\q-p^ l^Pj (4) 

where, || . || denotes the Euclidean norm. 

The Voronoi cell Vi is the collection of those points which 
are closest (with respect to the Euclidean metric) to pj 
compared to any other point in V . The boundary of each 
Voronoi partition is the union of a finite number of line 
segments forming a closed C° curve. In M^, the intersec- 
tion of any two Voronoi cells is either null, a line seg- 
ment, or a point. In d dimensional space, the bomidaries 
of the Voronoi cells are unions of convex subsets of at 
most d—1 dimensional hyperplanes in and the inter- 
section of two Voronoi cells is either a convex subset of 
a hyperplane or a null set. Each of the Voronoi cells is a 
topologically connected non-null set. 

Basic components of the Voronoi partition are 

i) A space which to be partitioned. 

ii) A set of sites, or nodes or generators. 

iii) A distance measure such as the Euclidean distance. 

3.2 Generalization of Voronoi partition 

Here we prcisent a generalization of the Voronoi partition 
considering the heterogeneity in the sensors' capabili- 
ties. We can find several extensions or generalizations of 
Voronoi partition to suit specific applications [9,18,19]. 
Herbert and Scidcl [24] have introduced an approach in 
which, instead of the site set, a finite set of real-valued 
functions : £) M is used to partition the domain D. 
Standard Voronoi partition and other known generaliza- 
tions can be extracted from this abstract general form. 

In this paper we define a generalization of the Voronoi 

partition to suit our application, namely the heteroge- 
neous multi-agent search. We use, 

i) The domain of interest as the space to be parti- 
tioned. 

ii) The site set as the set of points in the domain of 
interest which are the positions of the agents in it. 

iii) A set of node functions in place of a distance mea- 
sure. 

Consider a space Q C M'', a set of points called nodes or 

generators V = {pi,P2, ■ ■ ■ ,Pn}, Pi G Q, with pi pj, 
whenever i ^ j, and monotonically decreasing analytic 
functions fi : R"*" t—^ R, where fi is called a node func- 
tion for the i-th node. Define a collection {Vi}, i € 
{1,2,..., N}, with mutually disjoint interiors, such that 



Q = Ujg{i_2,...,jv}^, where Vi is defined as 

Vi = {q€Q\fi{\\pi-q\\)> fj{\\pj-q\\) ^ l, (5) 

je{i,2,...,7V}} 

We call {T^},iG{l,2,..., N}, as a generalized Voronoi 
partition of Q with nodes V and node functions fi. Note 
that 

( 1 ) Vi can be topologically non-connected and may con- 
tain other Voronoi cells. 

(2) In the context of the problem discussed in this pa- 
per, q &Vi means that the i-th agent /sensor is the 
most effective in sensing at point q. This is refiected 
in the > sign in the definition. In standard Voronoi 
partition used for the homogeneous case, < sign for 
distances ensured that i-th sensor is most effective 
inVi. 

(3) The condition that fi are analytic implies that for 
every i, j G {1, 2, . . . , N}, fi — fj is analytic. By the 
property of real analytic functions [25], the set of 
intersection points between any two node functions 
is a set of measure zero. This ensures that the inter- 
section of any two cells is a set of measure zero, that 
is, the boundary of a cell is made up of the union of 
at most d—1 dimensional subsets of R''. Otherwise 
the requirement that the cells should have mutu- 
ally disjoint interiors may be violated. Analyticity 
of the node functions fi is a sufficient condition to 
discount this possibility. 

Theorem 1 The generalized Voronoi partition depends 
at least continuously on V. 

Proof: If Vi and Vj are adjacent cells, then all the points 
q £ Q on the boundary common to them is given by 
{q e Q\ft{\\pi - q\\) = fjiWPj - g||)}, that is, the in- 
tersection of corresponding node functions. Let the j- 
th agent move by a small distance dp. This makes the 
common boundary between Vi and Vj move by a dis- 
tance, say dx. Now as the node functions are monoton- 
ically decreasing and are continuous, it is easy to see, 
that dx Q as dp ^ 0. This is true for any two i and 
j. Thus, the Voronoi partition depends continuously on 
'P = {pi,P2,...,Pn}- □ 



3.3 Special cases 

The name 'generalized Voronoi partition' suggests, that 
by suitably selecting parameters like the node functions, 
one should get the standard Voronoi partition and its 
generalizations as special cases. Below we discuss a few 
interesting special cases. 
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Case 1: Weighted Voronoi partition 



with Rp^ as a parameter specific to each node function. 



We consider multiplicatively and additively weighted 
Voronoi partitions as special cases. Let 



fiin) = -OiTi - di 



(6) 



where, ri =\\ Pi — q\\ and, ai and di take finite positive 
real values for i = 1, 2, . . . , A''. 

Thus, 

Vi = {q& QlaiTi + di < a^r^ + d^, Vj ^ i, (7) 
J e {1,2,. ..,7V}} 



The partition {Vi} is called a multiplicatively and ad- 
ditively weighted Voronoi partition, ctj arc called multi- 
plicative weights and di are called additive weights. With 
this generalization, the Voronoi colls no longer possess 
the nice property of being topologically connected sets. 
The Voronoi cells could be made up of disjoint patches 
and one or more Voronoi cells can get embedded inside 
another cell. 

Case 2: Standard Voronoi partition 



It can be noted that the Voronoi partitions with non- 
Euclidean metric or pseudo-metric, and with objects 
such as lines, curves, discs, polytopcs, etc., instead of 
points, as generators, can also be viewed as special cases 
of the generalized Voronoi partition (5). 



Case i 
tance 



Voronoi partition based on Non-Euclidean die- 



One of the ways in which a standard Voronoi partition 
can be generalized is use of non-Euclidean distance. The 
distance measure used depends on the application. As an 
example, in [15], authors use a non-Euclidian distance 
measure to incorporate anisotropy in sensor effective- 
ness. They use d{q,pi) = \\pi - q\\Li), with \\pi - = 
{q - PiVLi{q - Pi), where 



with 



Fi = 



Li = F/ Fi 



c/a \ / cos6'i sinOi 
c/b J y — sin^jcos^i 



(10) 



(11) 



The standard Voronoi partition can be obtained as a 
special case of (5) when the node functions are /i(r,) = 

Vi = {qG Q\\\pi - q\\ < \\p, - q\\ Vj / i, (8) 
jG{l,2,...,iV}} 



By setting the node function as fi{\\pi — in the 

generalized Voronoi partition (5), we get the generaliza- 
tion of Voronoi diagram used in [15]. 

Case 5: Other possible variations 



It can be shown that if the node functions are homoge- 
neous = /(.) for each i G {1, 2, . . . , N}), then the 
generalized Voronoi partition gives the standard Voronoi 
partition. 



Other possible variations of the Voronoi partition are 
using objects other than points as sites/nodes and gen- 
eralization of the space to be partitioned. It is easy to 
see that these generalizations can be obtained by specific 
site sets and the spaces. 



Case 3: Power diagram 

Power diagram PV or Voronoi diagram in Laguerre ge- 
ometry used in [13] is defined as 

PVi = {q€ Q\dp{q,pi) < dp{q,pj),i + j} (9) 

where, dp{q,pi) = \\pi — — R^., the power distance 
between q and pi , with Rp. > being a parameter fixed 
for a given node pi. In the context of robot coverage 
problem addressed in [13], Rp. represents the radius of 
foot print for i-th robot. It is easy to see that the power 
diagram can be obtained from the generalized Voronoi 
partition (5) by setting 

fi{\\q-Pi\\) = \\Pi-qf-Rl 



3.4 Generalized Delaunay graph 

Delaunay graph is the dual of Voronoi partition. Two 

nodes are said to be neighbors (connected by an edge) , if 
the corresponding Voronoi cells are adjacent. This con- 
cept can be extended to generalized Voronoi partition- 
ing scheme. For the sake of simplicity wc call such a 
graph a Delaunay graph, Qd. Note that the generalized 
Delaunay graph, in general, need not have the property 
of Delaunay triangulation, in fact, it need not even be a 
triangulation. 

Two nodes are said to be neighbors in a generalized, De- 
launay graph, if the corresponding generalized Voronoi 
cells are adjacent, that is, {i,j) € fg^, the edge set cor- 
responding to the graph Qoi'^^Vif^Vj 
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4 Heterogeneous locational optimization prob- 
lem 

Here we formulate and solve heterogeneous locational 
optimization problem (HLOP) for a mobile sensor net- 
work. Let (5 C ffi'* be a convex polytope, the space 
in which the sensors have to be deployed; ^ : Q i— > 
[0, 1] , be a continuous density distribution function; V = 
{pi,P2, ■ ■ ■ ,Pn}, Pi Q he configuration of N sensors; 
fi : R"*" R, i € {!,..., N}, be analytic, monotonically 
decreasing function corresponding to the «-th node, the 
sensor effectiveness function of i-th agent; and Vi G Q 
be the generalized Voronoi cell corresponding to the i-th 
node/sensor. 

The density (j){q) is the probability of an event of in- 
terest occurring in q € Q, indicating the importance of 
measurement at the given point in Q. As (f>(q) 1, the 
importance of measurement at q increases as the proba- 
bility of occurrance of an event of interest is higher. The 
objective of the problem is to deploy the sensors in Q so 
as to maximize the probability of detection of an event 
of interest. Let fi be the variation of the sensor effective- 
ness with the Euclidean distance. It is natural to assume 
fi to be monotonically decreasing. 

In case of homogeneous sensors, the sensor located in 
Voronoi cell Vi is closest to all the points q G Vi and 
hence, by the strictly decreasing variation of sensor's 
effectiveness with distance, the sensor is most effective 
within Vi. Thus, the Voronoi decomposition leads to op- 
timal partitioning of the space in the sense that, each 
sensor is most effective within its corresponding Voronoi 
cell. In the heterogeneous case too, it is easy to see 
that each sensor is most effective in its own generalized 
Voronoi cell, by the very definition of the generalized 
Voronoi decomposition. Now, as the partitioning is op- 
timal, we need to find the location of each sensor within 
its generalized Voronoi cell. 

4-1 The objective function 

Consider the objective function to be maximized, 
n{V)= [ m8^{fi{\\q-pi\\)}ct>{q)dQ 

Jq ' 

(12) 

i JVi 

where |[ . || is the Euclidean distance. Note that the gener- 
alized Voronoi decomposition splits the objective func- 
tion into a sum of contributions from each generalized 
Voronoi cell. Hence the optimization problem can be 
solved in a spatially distributed manner, that is, the op- 
timal configuration can be achieved, by each sensor solv- 
ing the part of objective function corresponding to its 
cell using only local information. 



Lemma 2 The gradient of the multi- center objective 
function (12) with respect to pi is given by 




where, ri = \\q-Pi\\. 

Proof. Let us rewrite (12) as 

n = Y,^i (14) 

i 

where, Hi = Jy, fi{ri)(j){q)dQ . Now, 

Applying the generalised form of the Leibniz theorem 
[26] 

/ 'i^i<l)?rin)dQ 
dpi Jvi dpi 

+ H / T^ij{q)-iiij{q)(l>{q)fi{ri)dQ (16) 




{q).Uji{q)(l){q)fj{rj)dQ 



where, 

(1) Ni is the set of indices of agents which are neighbors 
of the i-th agent in Gd, the generalized Delaunay 
graph. 

(2) Aij is the part of the bounding surface common to 
Vi and Vj. 

(3) iiij (g) is the unit outward normal to Aij at g € Aij . 
Note that nij{q) = —nji{q), \/q G Aij. 

(4) nij{q) = -^{q), the rate of movement of the 
boundary at g S Aij with respect to p,. Note that 
u^j{q) = Uji{q). 

(5) Note also that fi{ri) = fj{rj), Vg G Aij, by defini- 
tion of the generalized Voronoi partition. 

By (3)-(5) above, and as (p is continuous, it is clear 
that for each j e Ni, n^j{q).\iij{q)(j){q)fi{ri)dQ = 

- Iaj,, T^jM-^jii.<l)(t>{q)f3{rj)dQ 
Hence, 

□ 



6 



4-2 The critical points 

The gradient of the objective function (12) with respect 

to Pi, the location of the i-th node in Q, can be deter- 
mined using (13) (by Lemma 2) as 



OPi JVi "Pi 



m^ir^{Pi-(i)dQ 



Vi d{ri) 



-MvM-Cv.) 



(17) 



where, = and 4>{q) = -(f>{q)^^. As fi,i e 

{1, . . . , N} is strictly decreasing, (f>(q) is always non- 
negative. Here My. and Cy, arc interpreted as the mass 
and centroid of the cell Vi with (/) as density. Thus, the 
critical points are pi = Cy. , and such a configuration V, 
of agents is called a generalized centroidal Voronoi con- 
figuration. 

Theorem 3 The gradient, given by (1 7), is spatially dis- 
tributed over the Delaunay graph Qn- 

Proof. The gradient (17) with respect to pi £ P, the 
present configuration, depends only on the correspond- 
ing generalized Voronoi cell Vi and values of (f> and the 
gradient of fi within Vi. The Voronoi cell Vi depends 
only on the neighbors Mg^ [pi, P) oipt. Thus, the gradi- 
ent (17) can be computed with only local information, 
that is, the neighbors of Pi in ^/j. □ 

The critical points are not unique, as with the stan- 
dard Voronoi partition. But in the case of a generalized 
Voronoi partition, some of the cells could become null 
and such a condition can lead to local minima. 

4.3 The control law 

Let us consider the sensor dynamics as 

Pi = Ui (18) 

Consider the control law 

Ui = -kprop{Pi - Cvi) (19) 

Control law (19) makes the mobile sensors move toward 
for positive kprop- If, for some i G {1, . . . , N}, = 0, 
then we define Cy. = p,. 

It is not necessary that S Vi, but G Q is true 
always and this fact ensures that Q is an invariant set 
for (18) under (19). 

It is easy to see, that the control law (19) is spatially 
distributed in the generalized Delaunay graph. 



Theorem 4 The trajectories of the sensors governed by 
the control law (19), starting from any initial condition 
P{0), will asymptotically converge to the critical points 
ofH. 

Proof. Here we use LaSallc's invariance principle 
discussed earlier. Consider V{P) = —H , where 
P = {pi,p2, ...,Pn} represents the configuration of N 
agents/sensors. 



V{P) = - 



Itt 



^dpi 

: -2akprop ^ {pi -Cvif 



(20) 



We observe that, F : Q 1-^ K is continuously differen- 
tiable in Q as {V^} depends contnuously on P by The- 
orem 1; M = Q is a compact invariant set; V is nega- 
tive definite in M; = y-i(O) = {CyJ, which itself is 
the largest invariant subset of E by the control law (19). 
Thus by LaSalle's invariance principle, the trajectories 
of the agents governed by control law (19), starting from 
any initial configuration 7^(0), will asymptotically con- 
verge to set E, the critical points of H, that is, the gener- 
alized centroidal Voronoi partitions with respect to the 
density function as perceived by the sensors. □ 



5 Constraints on agents' speed 

We proposed a control law to guide the agents toward 
the critical points, that is, to their respective centroid, 
and observed that the closed loop system for agents mod- 
eled as first order dynamical system, is globally asymp- 
totically stable. Here we impose some constraints on the 
agent speeds and analyze the dynamics of closed loop 
system. 

5. 1 Maximum speed constraint 

Let the agents have a constraint on maximum speed of 
Umaxii for i = 1, . . . , n. Now consider the control law 

I kpropiPi ^ C'v'j), If \kprop{Pi C'y^ ) | < Umaxi 
Ui = \ iPi-Cvi) 



-Urr 



Ml(pi-Cn) 



Otherwise 



(21) 



The control law (21) makes the agents move toward their 
respective centroids with saturation on speed. 

It is easy to see, that the control law (21) is spatially 
distributed in the generalized Delaunay graph. 
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Theorem 5 The trajectories of the agents governed by 
the control law ( 21 ), starting from any initial condition 
'P(O) e Q^, will asymptotically converge to the critical 
points ofH^. 

Proof. Consider = — W, whereP = {p\,p2, ■■■,Pn} 

represents the configuration of N agents. 



V{V) 



dH 
~dt 



E 



dpi 



-Pi 



(22) 



ie{l,2,...,N} 

'-2aX]iG{l,2,...,Af} MviiPi - Cv,)kprop 
X {Pi - Cvi), If \kprop{Pi -Cvi)\<U„ 

-2aEie{i,2,...,Jv}^v'. 

X {jpi - CVi)Umaxt - ^ 

otherwise 



''(\\P^-CV^ 



-2akprop Eiefi,2,...,iV} ^vMPi - C'v, ||)^ 

If \kpj,op{Pi < Umaxi 

— 2a Eie{l,2,...,Ar} ^maxi^Vi 



\\Pi-Cv 



otherwise 



l|Pi-CvJ 



We observe that F : Q i— > M is continuously differen- 
tiable in Q as {Vi} depends at least continuously on V 
(Theorem 1), and V is continuous as u is continuous; 
M = Q is a compact invariant set; V is negative def- 
inite in M; E = ^"^(0) = {CyJ, which itsch is the 
largest invariant subset of E by the control law (21). 
Thus, by LaSalle's invariance principle, the trajectories 
of the agents governed by control law (21), starting from 
any initial configuration 'P(O) G Q^, will asymptotically 
converge to the set E, the critical points of H, that is, 
the generalized ccntroidal Voronoi configuration with re- 
spect to the density function as perceived by the sensors. 

□ 



5.2 Constant speed control 

The agents may have a constraint of moving with a con- 
stant speed Ui . But we let the agents slow down as they 
approach the critical points. For i = 1, . . . , n, consider 
the control law 



if \\Pr - CvA\ > S 



Ui{pi — Cvi)/5, otherwise 



(23) 



where, 5 > 0, predefined value, such that the control 
law (23) makes the agents move toward their respective 
centroids with a constant speed of Ui when they at a 
distance greater than 5 from the corresponding centroids 
and slow down as they approach them. 



It is easy to sec, that the control law (23) is spatially 
distributed in the generalized Delaunay graph. 

Theorem 6 The trajectories of the agents governed by 
the control law (23), starting from any initial condition 
V{0) G , will asymptotically converge to the critical 
points ofH. 

Proof. Consider V{V) = —H, where V = {pi,P2, ■■■,Pn} 
represents the configuration of N agents. 



V{V) = - 



Itt 



E 



ie{l,2,...,JV} 



dpi 



-Pi 



(24) 



if \\Pi - CvA\ > s 

-2aE.e{i,2,...,iv} UiMvM - Cv^/d, 
otherwise 



We observe that y : Q i-^ K is continuously differen- 
tiablc in Q as {Vi} depends at least continuously on V 
(Theorem 1), and V is continuous as u is continuous; 
M = Q is a compact invariant set; V is negative def- 
inite in M; ^ = y-i(O) = {CyJ, which itself is the 
largest invariant subset of E by the control law (23). 
Thus, by LaSalle's invariance principle, the trajectories 
of the agents governed by control law (23), starting from 
any initial configuration ■p(O) € Q^, will asymptotically 
converge to the set E, the critical points of 7i, that is, 
the generalized centroidal Voronoi configuration with re- 
spect to the density function as perceived by the sensors. 

□ 



6 Limited range sensors 

In reality the sensors will have limited range. In this sec- 
tion we formulate a spatially distributed limited range 
locational optimization problem. 

Let Ri be the limit on range of the sensors and B{pi, Ri) 
be a closed ball centered at pi with a radius of Ri. The 
i-th sensor has access to information only from points in 
the set ViT] B{pi, Ri). Let us also assume that fi{Ri) = 
fj{Rj), Vi,j e {!,..., N}, that is, we assume that the 
cutoff range for all sensors is the same. Consider the 
objective function to be maximized, 

nV) = ^ / MQ)k\\Pi - Q\\))dQ (25) 

^ J{VinB(pi,R)) 

'f,{r) if r<Ri 



where, fi{r) 



fi{Ri) otherwise 
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2 4 6 8 10 

Fig. 1. Illustration of fi and fi in the presence of a limit on 
sensor range. The solid curve represent the sensor effective 
function fi and dotted curve is the actual sensor effective 
function /,(r) = fi{r) + [1 - fi{R)) with i? = 6. 

It can be noted that the objective function is made up 
of sums of the contributions from sets Vi n B{pi,Ri), 
enabling the sensors to solve the optimization problem 
in a spatially distributed manner. 

In reality for range limited sensors the effectiveness 
should be zero beyond the range limit. Consider 

/i(-) = fi{-) ~ Ii{Ri) (Figure 1). It can be shown that 
the objective function (25) has the same critical points 
if fi is replaced with ft, as the difference in two objec- 
tive functions will be a constant term fi{Ri) (Note that 
we have assumed fi{Ri) = fj{Rj)Vi,j e {1, . . . , N}.). 

The gradient of (25) with respect to pi can be determined 
as 

— — (P) = 2M(y^n^(p. j^))(C(y.p,g(p. jj) - Pi) (26) 



We use the control law 

Ui = -kpropiPi - C'(VinB(p.,fl))) (27) 



It is easy to show, that the gradient (26) and the control 
law (27) are spatially distributed in the r-limited Delau- 
nay graph Qld, the Delaunay graph incorporating the 
sensor range limits. 

Theorem 7 The trajectories of the sensors governed by 
the control law (27), starting from any initial condition 
■p(O), will asymptotically converge to the critical points 
ofU. 

Proof. The proof is similar to that of the Theorem 4 with 
V = —H{V). It is It can be shown that V is continu- 
ously differentiable based on Theorem 2.2 in [12], where 
a similar result is proved. □ 



2 4 6 8 10 2 4 8 8 

(a) (b) 
Fig. 2. Trajectories of sensors along with final Voronoi 
cells, with (a) an uniform uncertainty density, and (b) 

(l,{x,y) = 0.9e'' °*«"-"''+(''-i»''). indicate initial posi- 
tions of the sensors, and 'o' indicate the centroid of Voronoi 
cells corresponding to agent locations at the end of simula- 
tion. 

7 Illustrative example 

Here we provide simulation results as an illustration 
of the optimal deployment using the proposed het- 
erogeneous locational optimization. In this example 
we consider 10 mobile sensors without any sensor 
range limits in a square area of size 10 x 10. We used 
fi = —ctirf as node functions, with {ai, . . . aio} = 
{1, 1.25, 1.5, 0.75, 0.8, 1.3, 0.9, 1.1, 1.4, 1.2}. Figure 2 (a) 
shows the trajectories of the sensors with a uniform un- 
certainty density and Figure 2 (b) shows the trajectories 
with (j){x,y) = Q^<^e-°m(x-wf+{y-w?)^ The sensors 
reach the optimal configuration which is the centroids 
of multiplicatively weighted Voronoi cells in this case. 
In both figures, 'o' represents the centroids of the gen- 
eralized Voronoi cell corresponding to the final agent 
configuration, and it can be seen that the each agent is 
sufficiently close to corresponding centroid. The simu- 
lations were terminated when max^ \\pi — (C)\/. || < 0.5. 
The generalized Voronoi partition corresponding to 
final agent configurations are also shown in Figure 2. 

Figure 3 (a) shows the plot of the normalized objective 
function with time steps. It can be observed that the 
control law successfully makes the agents/sensors move 
toward the optimal deployment maximizing the objec- 
tive function. As noted earlier, the strategy does not 
guarantee global optimal configuration. Plot of an error 
measure of 2 N^Pi ~ {^)viY with time steps 

is shown in Figure 3 (b). It can be observed that the 
control law successfully drives the agents toward the re- 
spective centroids. 

8 Conclusions 

A generalization of Voronoi partition has been pro- 
posed and the standard Voronoi decomposition and its 
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5 10 15 20 25 30 5 10 15 20 25 30 

Time Steps TimeSleps 



(a) (b) 

Fig. 3. The history of (a) normahzed objective function and 
(b) normalized error measure 

variations are shown as special cases. A heterogeneous 
spatiaUy distributed locational optimization has been 
formulated and solved using the proposed partitioning 
scheme and the a variation of centroidal Voronoi con- 
figuration is shown to be the optimal deployment of the 
limited range sensors with heterogeneous capabilities. 
Illustrative simulation result was provided to support 
the theoretical results. 

References 

[1] C. G. Cassandras and W. Li, Sensor networks and cooperative 
control, European J. Control, vol. 11, no. 4-5, 2005 
pp. 436-463. 

[2] W. Li and C.G. Cassandras, Distributive cooperative 
coverage control of mobile sensing networks, Proc. IEEE 
Conference on Decision and Control, and the European 
Control Conference 2005, Seville, Spain, December, 2005, 
pp. 2542-2547. 

[3] Y. Zou and K. Chakrabarty, Sensor deployment and 
target localization baesd on virtual forces, Proc. of IEEE 
INFOCOM, 2003, pp. 1293-1303. 

[4] LI. Hussein and D.M. Stipanovic, Effective coverage control 
using dynamic sensor networks with flocking and guaranteed 
collision avoidance. Proceedings of the 2007 American Control 
Conference, New York, USA, July, 2007, pp. 3420-3425. 

[5] LI. Hussein and D.M. Stipanovic, Effective coverage control 
using dynamic sensor networks. Proceedings of the 45th IEEE 
Conference on Decision and Control, San Diego, CA, USA, 
December, 2006, pp. 2747-2752. 

[6] I. I. Hussein and D. M. Stipanovic, Effective coverage 
control for mobile sensor networks with guaranteed 
collision avoidance, IEEE Transactions on Control Systems 
Technology, vol. 15, no. 4, July 2007, pp. 642-657. 

[7] I.I. Hussein, A Kalman filter based control strategy 
for dynamic coverage control, Proc. of the American 
Control Conference (ACC 2007), New York, July, 2007, 
pp. 3271-3276 

[8] Z. Drezner, Facility Location: A Survey of Applications And 
Methods, New York, NY: Springer, 1995. 

[9] A. Okabe and A. Suzuki, Locational optimization problems 
solved through Voronoi diagrams, Eurpean Journal of 
Operations Research, vol. 98, no. 3, 1997, pp. 445-456. 



[10] Q. Du, V. Faber and M. Gunzburger, Centroidal Voronoi 
tessellations: applications and algorithms, Siam Review, Vol 
41, no. 4, 1999, pp. 637-676. 

[11] Guruprasad K. R. and D. Ghose, Deploy and search strategy 
for multi-agent systems using Voronoi partitions. Proceedings 
of Jfth international Symposium on Voronoi Diagrams 
in Science and Engineering (ISVD 2007), University of 
Glamorgan, Wales, UK, July 9-11, 2007, pp. 91-100. 

[12] J. Cortes, S. Martinez, and F. BuUo, Spatially-distributed 
coverage optimization and control with limited-range 
interactions, ESAIM: Control, Optimization and Calculus of 
Variations, vol 11, no. 4, 2005, pp. 691-719. 

[13] L.C.A. Pimenta, V. Kumar, R.C. Mesquita, and A.S. Pereira, 
Sensing and coverage for a network of heterogeneous robots, 
to appear in Proc. of IEEE Conference on Decision and 
Control, 2008. 

[14] K.R,. Guruprasad and D. Ghose, Heterogeneous sensor 
based Voronoi decomposition for spatially distributed limited 
range locational optimization, in Voronoi's Impact on 
Modern Science, Book 4t 2, Proceedings of 5th 

Annual International Symposium On Voronoi Diagrams 
(ISVD 2008) Eds., D.-S. Kim, K. Sugihara, Kiev, Ukraine, 
September, 2008, pp. 78-87. 

[15] A. Gusrialdi, S. Hirche, T. Hatanaka, and M. Fujita, Voronoi 
based coverage control with anisotropic sensors. Proceedings 
the American Control Conference, 2008, pp. 736-741. 

[16] G. Voronoi. Nouvelles applications des paramtres continus 
la thorie des formes quadratiques. Journal fr die Reine und 
Angewandte Mathematik, vol. 133, 1907. pp. 97-178. 

[17] G.L. Dirichlet. ber die Reduktion der positiven quadratischen 
Formen mit drei unbestimmten ganzen Zahlen. Journal fr 
die Reine und Angewandte Mathematik, vol. 40, 1850, pp. 
209-227. 

[18] F. Aurenhammer, Voronoi diagrams - A survey of a 
fundamental geometric gata structure, ACM Computing 
Surveys, vol 23, no. 3, 1991, pp. 345-405. 

[19] P.A. Arbelaez and L.D. Cohen, Generalized 

Voronoi tessellations for vector-valued image segmentation, 
Procedings of 2nd IEEE Workshop on Variational, Geometric 
and Level Set Methods in Computer Vision (VLSM'03), 
September 2003, Nice, France. 

[20] R. Diestel, Graph Theory 2nd edn., Ser. Graduate Texts in 
Mathematics, New York, Springer- Ver lag 173 (2000). 

[21] J. P. LaSalle, Some extensions of Liapunov's second method, 
IRE Transactions on Circuit Theory, CT- vol. 7, no. 4, 
December, 1960, pp. 520-527. 

[22] E.A. Barbashin and N.N. Krasovski, Ob ustoichivosti 
dvizheniya v tzelom, Dokl. Akad. Nauk., USSR, vol. 86, no. 
3, 1952, pp. 453-456. 

[23] H.J. Marquez, Nonlinear Control Systems - Analysis and 
Design, John Wiley & Sons, Inc., 2003. 

[24] H. Edelsbrunner and R. Seidel, Voronoi diagrams and 
arrangements. Discrete Comput. Geom., vol 1, 1986, 
pp. 25-44. 

[25] S.G. Krantz and H.R. Parks, A Primer of Real Analytic 
Functions (2nd edn.), Birkhauscr Advanced Texts, Basler 
Lehbiicher, 2002. 

[26] P.K. Kundu and I.M. Cohen, Fluid Mechanics (2nd edition). 
Academic Press, 2002. 



10 



